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ABSTRACT 

We consider theoretical models of emission of TeV photons by Cyg X- 1 during a flare discov- 
ered by the MAGIC detector. We study acceleration of electrons to energies sufficient for TeV 
emission, and find the emission site is allowed to be close to the black hole. We then consider 
pair absorption in the photon field of the central X-ray source and a surrounding accretion 
disc, and find its optical depth is < 1, allowing emission close to the black hole. On the other 
hand, the optical depth in the stellar field is ~10 at ~1 TeV. However, the optical depth drops 
with increasing energy, allowing a model with the initial energy of > 3 TeV, in which photons 
travel far away from the star, initiating a spatially extended pair cascade. This qualitatively 
explains the observed TeV spectrum, though still not its exact shape. 

Key words: accretion, accretion discs - radiation mechanisms: non-thermal - stars: individ- 
ual: Cyg X-l - stars: individual: HDE 226868 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

The black-hole binary Cyg X-l is at present the only X-ray binary 
that is both certainly accretion-powered and from which very high 
energy y-rays, >0.1 TeV, have been detected by one of the cur- 
rent generation of TeV detectors. Nam ely, a transient TeV emission 
was detected by the MAGIC telescope dAlbert et alj2007l . hereafter 
A07). Among the three other currently known TeV-emitting bina- 
ries (which are all per sistent in the TeVs ex cept for orbital varia- 
tions), PSR B 1259-63 jAharonian et al.l2005h is known to be pow- 
ered by the rotation energy of a young pulsar rather than by ac- 
cretion jjohnston et al.llT992h . with the resulting interaction of the 
pulsa r and stellar winds leadin g to production of en ergetic pho- 
tons dMaraschi & Treveslll98lh . Then, LS I +61° 303 ^Albert et all 
l2006h is very likely of the same nature, as indicated by its radio 
observation s showing a cometar y-tail-like variable structure rather 
than a jet dPhawan et alj [2006). A possible further evidence for 
the presence of a young pulsar in that system is provided by a 
detection of a powerful ~0.2-s flare with the peak 15-150 keV 
flux of =! 5 x 10~ 8 erg cnr 2 s~ l (higher by a factor of ~ 10 3 
than the aver age source flux) and a blackbody-like spectrum by the 
SWIFT BAT dPe Pasquale et alj2008llBarthelmv et alj2008h . Such 
an event is consist ent with a soft gamma-r ay repeater/anomalous 
X-ray pulsar burst JDubus & Gie bels 2008). Then, the case of LS 
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5039 jAharonian e t al. 2006) is least clear, but the striking similar- 
ity of the broad-band spectra of all three sys tems points to similar 
mechanisms of their activity (Dubus 2006b). All three persistent 
systems are also similar in being eccentric high-mass binaries con- 
taining massive B or O stars. Furthermore, their broad-band spectra 
are very much unlike those of accreting sources at similar Edding- 
ton ratios as well as are unlikely to b e dominated by a col limated 
jet emission (see, e.g., a discussion in IZdziarski et al j2008h . 

Then, there were some early claims of TeV emission from ac- 
creting X-ray binaries, most notably from Cyg X-3. Its TeV flux 
was claimed to be very strong and modulated at the orbital period 
(e.g.. lLambetaflll982h . which, however, was not confir med by a 
subsequent Whipple Observatory pointing jWeekes| [l983l). Cyg X- 
3 has not, until now, been detected by either the MAGIC, VERI- 
TAS, or HESS telescopes. 

Thus, the case of TeV emission from Cyg X-l is of crucial 
importance for our understanding of emission of y-rays from ac- 
cretion flows. The statistical significance of the detected transient 
was 4.9<x or 4. lcr post-trial (A07), unlikely to appear by chance. 
Supporting the reality of this transient emission is als o the coinci- 
dence with a flare of Cy g X-l seen in X-rays (A07: iTurler et al.l 
l2006l:lMalzac et al.ll2008l , hereafter M08). 

The TeV detection occurred on 2006 September 24, MJD 
54002. [928-982] (A07), which corresponds to the orbital phases 
of 0.90-0.91, i.e., with the black hole behind the companion. No 
signal was detected (A07) during the first half of that MAGIC ob- 
servation, MID 54002.[875-928], indicating time variability of the 
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Figure 1. The spectral components of the Cyg X-l/HDE 226868 system. 
The long dashes give the blackbody spectrum of the OB supergiant (for 
the parameters of Herrer o et alj l 1995i). The dotted curve gives the thermal 
Comptonization component from the hot inner flow. The short dashes give 
the spectrum of the surrounding optically-thick disc (tidally truncated at the 
outer radius), and the solid curve gives the spectrum of the disc together 
with the hot-plasma emission. See Section[2]for details. 



TeV emission on the time scale of less than an hour. The emission 
was detected in the 0.15-1 TeV range, and was fitted by a power 
law with a photon spectral index of T = 3.2 ± 0.6 (A07). The 0.1-1 
TeV isotropic luminosity corresponding to the be st fit is ^ 2 x 10 34 
erg s -1 (assuming the distance of D = 2 kpc, see Zi61 kowski1l2005l 
and references therein). The MAGIC telescope also observed, but 
did not detect, Cyg X-l on 25 other nights, and obtained an upper 
limit to the flux vs. energy an order of magnitude below the flux 
during the flare (A07). 

As mentioned above, the TeV flare occurred during a longer 
X-ray outburst. Based on data from the RXTE All Sky Monitor, the 
X-ray flare lasted three days, about MJD 54002.0-54005.0. Cyg X- 
1 was then in the hard spectral state (M08). The INTEGRAL IBIS 
and SPI instruments observed Cyg X-l in the 18-700 keV range 
on MJD 54002.4-54004.3 (M08), during which the estimated bolo- 
metric luminosity was a 4 x 10 37 erg s -1 . M08 obtained, in partic- 
ular, the X-ray spectrum for an interval coinciding with the TeV 
detection. The spatial location of the TeV emission region remains 
unknown. However, the simplest hypothesis to consider that it was 
close to the X-ray source. Then, we can calculate e* pair absorption 
of the TeV photons as well as the Compton energy loss of acceler- 
ated electrons due to interactions with the X-rays. 



2 THE BINARY PARAMETERS AND THE 1 eV-1 MeV 
SPECTRUM 

The system parameters of interest to our study are the stellar tem- 
perature, r*, the stellar radius, R+, the separation, a, and the in- 
clination, i. The radius of the companion, HDE 226868, and the 
compone nt masses are subject to a considerable uncerta inty, com- 
pare, e.g., iHerrero et all dl995h with Ziolkowskil J2005h . Here, we 
adopt the best model of lHerrero et al.l Jl995h . with R* = 1.2 x 10 12 
cm, the separation of a = 2A1R*, and T* ^ 3.2 x 10 4 K. The ob- 
served stellar flux is plotted in Fig. [T] The value of the inclination 
of Cyg X-l is relati vely uncertain; here, we adopt i = 30° (e.g., 
lGies&Boltonlll986t) . 



The 18-700 keV X-ray spectrum during the TeV flare was 
fitted by M08 as an e-folded power law. Physically, it corresponds 
to thermal Comptonization of blackbody photons produced at the 
inner edge of a cold disc. At low energies, the disc spectrum has a 
low-energy break, with the Rayleigh- Jeans form below the break. 
To model it, we assume the observed broad-band energy flux, F h , 
is given by, 



EF h (E) 



K(E b /lkeV) 2 - r exp(-E/E c ) 
(£/£ b )- 3 + (£/£ b ) r - 2 ' 



(1) 



which has the form of KE exp(- E/E c ) at E » E^. We assume 
Eb = 0.3 keV dFrontera et all2001D . For the remaining parameters, 
we use the best-fit parameters of M08, T = 1.34, E c = 111 keV, 
K = 1.25 keV crrT 2 s~ x . This spectrum is plotted in Fig.Q] 

The corresponding photon density within the central hot 
source, n h (E), can be calculated by assuming an isotropic and ho- 
mogeneous spherical source with the radius, R^, in which the aver- 
age photon escape time is 3^h/4c, 



E z n h (E,R h ) 



9D 2 EF h (E) 
ARTc 

a 



(2) 



(3) 



The hard state emission is well modelled as coming from a hot inner 
flow surrounded by a cold dis c truncated at ~30 -100^,,, where R g is 
the gravitational radius (e.g.. |Pone et al ] l2007h . For the black hole 
mass of M x ~ IOMq, the radius of the inner hot flow is R h ~ 10 s 
cm, which value we assume hereafter. 

Moving away from the hot plasma, the local photon den- 
sity includes an increasing number of soft photons from the sur- 
rounding disc. For a given height above the disc, most of the pho- 
tons come from the part of the disc within the radius equal to 
the height. The inner disc temperature is approximately given by 
2£7bb,max ~ Eb — 0.3 keV, with which factor the spectrum is expo- 
nentially cut off at high energies. The temperature decreases with 
the increasing radius as T hb oc R~ 3/4 , and the resulting disc black- 
body spectrum is Fa(E) oc E ,/3 . At the outer edge of a considered 
region, R out , the spectrum changes to the Rayleigh- Jeans one. The 
corresponding energy is E out = (R out /R h y 3l4 E h . Thus, we approxi- 
mate the disc spectrum by, 

^ 3 exp(-£/E b ) 
d( ' ° u,) " (E/E out y 3 + (E/E mt )-*' 3 ' 

where the normalization constant, K' , follows from the condition 
that Fd and F^ intersect at some energy, which we choose as 1 keV 
based on the fits to BeppoSAX data in Frontera et al. (2001). This 
choice also reproduces the soft excess appearing in the hard state 
of Cyg X-l at E < 3 keV dGierlinski etal]|l997l) . We assume the 
disc is tidally truncated at 10 12 cm. The spectrum of equation {3} 
for R out = 10' 2 cm and the total spectrum of the disc and the hot 
flow are plotted in Fig. [T] The disc photon density, n&(E, R m), is 
approximated analogously to equation 

E 2 n d (E, R mt ) - 9D 2 £F d (£, J R 0Ut )/(4^ 2 t f). (4) 



3 ELECTRON ACCELERATION 

Here, we consider whether electrons can be accelerated to TeV en- 
ergies within or close to the central X-ray source. We consider the 
acceleration rate in a generic form, applicable, e.g., to shock accel- 
eration, 



7ac 



1.8 x 10 7 £— s- 
1 o 



(5) 
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Figure 2. The Compton energy loss rate, — y, in the X-ray photon field, 
given by equations (TJ2J and assuming the radius of 10 9 cm (solid curve), 
and in the stellar field diluted at the position of the black hole (long dashes). 
The synchrotron loss rate is shown for B = 20 G (short dashes), and for 
B = 5 G (dots), and the dot-dashed line is (-) the acceleration rate, equation 
E) for B = 20 G and £ = 0.5. 



where y is the Lorentz factor, R^ = ym e c 2 /eB =; 1700y/B cm is 
the Larmor radius, B is the magnetic field strength, e is the electron 
charge, and £ < 1 is a dimensionless quantity parameterizing the 
acceleration efficiency (see, e.g., Bednarek & Giovannellfo oj 7 !). 
The maximum electron Lorentz factor can be then obtained for the 
acceleration rate being equal to an energy loss rate. 
The synchrotron energy loss is given by, 



, , , /(-£)' S : - (6) 

3 ra c c 8i \1G/ 

Equating this loss rate to the acceleration rate, equation (5), we 
obtain the corresponding maximum Lorentz factor of 

" 2 8 1/2 / B Y 112 



-7s 



/ max,S 



{o-tB) 



1.2 x 10 8 



Note that this implies a universal value of th e maximum energy of 
the synchrotr on component, as po inted out bv lGuilbert et al.l ( 1983) 
(and later bv lde Jager et alj|l996t) . Taking the characteristic syn- 
chrotron energy of (3/4n)y 2 (heB/m c c), we obtain, 



= - 40fMeV 



(8) 



16 at 

(where a f is the fine-structure constant), which is well below the 
TeV range. 

In addition, if electrons are accelerated within the hot plasma 
region, the Larmor radius has to be smaller than the acceleration 
site size, /? acc , which yields the maximum Lorentz factor of 

BRxx-e , ,„*B R. c 



Tmax.L 



6 x 10 3 



(9) 



ot c c- lG10 9 cm' 

Then, comparing equations Q and (9), we obtain the optimum 
magnetic field yielding the highest maximum the Lorentz factor 
subject to both constraints together, 



(10) 



Note that there is a spurious factor of n in eq. (3) of 

iBednar ek & Giovannelli (2003), an£ i me — 2 factor in their eq. (6) 
should be outside the argument of the logarithm. 



which Lorentz factor equals 

1/3 

=! 2 x 10 7 £ 1/3 



y, 



6ne 2 ^R x , 
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ra c c z cr T 



1 10" cm) ' 
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We see that the optimum magnetic field is relatively weak, e.g., 
= 160£ 1/3 G if Bacc = Rh = 10 s cm. The magnetic field expected 
in the hot inner plasma from equipartition arguments is consider- 
ably higher, B> 10 5 G (e.g.. iMalzac & Belmontll2008l) . Thus, the 
preferred acceleration site is outside of the hot inner plasma. It can 
be, e.g., a shock region where the stellar wind collides with the 
hot inner flow, or an inner part of the jet. The Lorentz factor of 
the accelerated electrons is completely sufficient to account for the 
photon energies observed in the TeV flare. In fact, the maximum 
Lorentz factor can still be lower than that of equation i ll It , which 
will then allow a value of the magnetic field higher than that of 
equation d 1 0b . 

An important process is also the loss rate due to Compton 
scattering with the X-ray and disc photons, equations lQJ|4j, and 
due to the stellar photon density at the position of the black hole. 
For simplicity, we assume isotropic photon fields. The energy loss 
of electrons with y 2 » 1 in s cattering with isotropic photons of 
an arbitrary ene rgy is given by I Jonesl d 19681) (see also eq. A30 in 
Zd ziarski|[l988T) . We integrate this rate over the photon density at 
the radius of 10 9 cm, obtaining the —y shown in Fig. [2] which cur- 
vature is due to to Klein-Nishina effects. We see that if the magnetic 
field is given by equation l |10t , the Compton energy losses rate are 
less than those due to the synchrotron losses. Thus, they would only 
slightly affect our estimates above. Also, we see that the losses on 
the X-rays dominate over those on the stellar photons. 

Thus, we find the conditions around the X-ray source to al- 
low electron acceleration up to y ~ 10 7 , which then easily allows 
emission of photons in the TeV range. The acceleration time scale, 
y/7acc> is then ver y short, < 1 s, see equation 



4 PAIR ABSORPTION OF TeV EMISSION 

Pair absorption of the TeV photons in photon-photon e* pair pro- 
duction events is a major issue in interpretation of the TeV flare. 
In a study of pair absorp t ion by photons from the companion star, 
IBednarek & Giovannellil d2007l) found that TeV photons should be 
preferentially observed from Cyg X-l at the orbital phases of 0.45- 
0.55, i.e., near the inferior conjunction, which is almost exactly op- 
posite to the phase of the observed event. However, pair absorption 
can also take place on size scales much smaller than the system sep- 
aration, due to pair-production events with the intense X-ray field of 
the central source. Below, we consider separately interaction with 
the X-rays and with the stellar photons. 



4.1 TeV-X-ray pair absorption 

As discussed above, the acceleration site is most likely larger than 
the hot inner plasma. Still, for completeness we calculate the e* 
pair absorption optical depth across the inner hot plasma, with the 
size of R\^, 



T h (E y ,R h ) 



Rh f n 



n h (E)a- y J(E y E) ll2 ]dE oc R 



(12) 



where o" vv is the cross section averaged ov er isotropic target p ho- 



;t p h 

tons (Gould & Schred e3l967l; correct ed in lBrown et al.l 1973bl : see 
also eqs. B11-B13 in Zdziarski 1988). The resulting optical depth 
is shown in Fig. [3] We see that whereas it is > 1 in the ~0.01-100 
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Figure 3. The optical depth to photon-photon pair production in the field 
of the hot plasma at R = 10 8 cm (dolid curve), and of the disc emission. 
The dotted and dashed curves give the optical depth due to the disc photons 
along the disc axis from the height of 10 8 and 10' cm, respectively. The 
assumed inner disc radius is 10 8 cm and the inner temperature is kT^ m3[X = 
0.15 keV. The acceleration site is likely at > 10' cm, Section[5] and thus the 
TeV photons are only weakly attenuated in the field of the disc and the hot 
plasma. 

GeV range, the (Klein-Nishina) high-energy decline of the cross 
section results in T h < 1 above this range. 

We then conside r pair absorption due to the disc photons. We 
follow the method of Bednarek (1993). The disc blackbody spec- 
trum is taken as having the temperature of fcr bbmax = 0.15 keV at 
the disc inner edge of R h = 10 8 cm, and diluted (due to a colour cor- 
rection; the normalization is also affected by the exact value of R^, 
which is unknown) to match the level of the disc spectrum shown 
in Fig. [T] We calculate the optical depth along the disc axis from 
the height of 10 s and 10 9 cm, see Fig. [3] We see that even in the 
former case, it is < 1 at E > 0. 1 TeV. For the distance of 10 9 cm 
(see Section[3}, pair absorption of photons > 0. 1 TeV on either the 
disc or hot-plasma photons is negligible. 

4.2 Pair absorption on stellar photons and spatially extended 
pair cascades 

We then calculate the optical depth on stellar photons, taking into 
account the ir anisotropy out side the stellar surface. We use the 
method s of iBednarekl i 19971) (with triple integration) and Dubus 
(2006a) (with quadruple integration) and found they fully agree 
with each other. The results for the orbital phase of 0.9, the bi- 
nary parameters given in SectionfZ] and three values of ( are shown 
in Fig.g] We see that at i = 30°, t, > 1 between 40 GeV and 20 
TeV. This might have implied that the acceleration process could 
not have taken place in the black hole vicinity. 

However, we note that a similar problem is faced in the case 
of the TeV-emitting X-ray binary LS 5039. Its parameters, R, m 
6 x 10" cm, T, 3.9 x 10 4 K, a =* 2.2R,, and the eccentricity of 
e =- 0.35 dCasares et af]|2005h . imply a system even more compact 
and, especially at periastron, more opaque to pair absorption than 
Cyg X-l. Still, a steady TeV emission is detected from LS 5039 at 
both the superior conjunction and the periastron dAharonian et al.l 
2006), which two points occur very close from each other in the 
orbital phase in that system. This is contrary to the theoretical pre- 
dictions based on the pair absorption alo ne, giving fluxes at those 
phase s much lower than those observed jDubusll2006al ; IBednarekl 
2006), similarly to the predictions implied by Fig.[4] The solution to 
this dilemma may be presented by pair cascades. They may be ini- 
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Figure 4. The pair-production optical depth in the field of the star from the 
location of the black hole up to infinity at the orbital phase of 0.9 and for 
the inclinations of 30°, 45° and 60°, from bottom to top. 
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Figure 5. The INTEGRAL spectrum from the ISGRI (cyan crosses), PICsIT 
(the cyan cross at the highest energy) and the SPI (magenta crosses), from 
M08, simultaneous with the MAGIC spectrum (black crosses, A07). The 
MAGIC spectrum is compared to the results from our pair cascade models, 
with the green and histograms corresponding to the primary monoenergetic 
electron injection at 3.16 TeV and 10 TeV, respectively. The blue histogram 
is for 3.16 TeV, but with the stellar temperature assumed to be 3.0 X 10 4 K 
(which is within the observational uncertainties for Cyg X-l). The normal- 
ization of the histograms have been chosen to match the MAGIC data. 

tiated by electrons and photons at energies of ~ 3- 10 TeV, at which 
energy the pair optical depth measured from the compact object 
drops to ~1, implying that those photons can propagate relatively 
far from the star. The resulting pairs produce photons at lower ener- 
gies, which then can propagate to the observer in the much diluted 
blackbody field. Such spatially extende d pair cascade model was 
applied to LS 5039 by Bednarek (20071) . Similar pair cascades can 
also take place in the X-ray field (see above). 

To calc u late p air cascades, we use two methods. One is 
of Bednarek] jl997l) , which uses an approximation in the Klein- 
Nishina regime, and the other is an independent Monte Carlo 
method using the full Klein-Nishina cross section. We have found 
the two methods give results consistent with each other. We present 
some of our results of the latter in Fig. [5] for monoenergetic elec- 
tron injection at 3.16 TeV and 10 TeV. This form of injection is 
supported by re sults of studies o f acceleration including energy 
losses, see, e.g., IProtherod Jioolh . We find we can reproduce the 
MAGIC spectrum only qualitatively. Generally, our models predict 
too much emission at energies > 1 TeV, and too little at < 0.2 TeV. 
We have found this to be a generic problem, and changing the in- 
jection energy does not improve the agreement with the data. Some 
improvement at the low energy end appears when a steep power 
law injection is added; however, this is at the expense of greatly in- 
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creasing the intrinsic power of the flare. On the other hand, we con- 
sider it possible that, given complications in interpreting Cherenkov 
telescope data, the actual spectrum of Cyg X-l was similar to that 
predicted by our models. At < 10 GeV, the pair cascade process 
stops, and the spectrum has the T = 1.5 photon index, as expected 
for an injection at high energies. The 0.1-1 TeV isotropic power is 
0.011-0.018 of the the total power injected in TeV electrons, for 
the three models shown in Fig.[5] This implies an injected power of 
~ 10 36 erg s -1 , which is still much less than the bolometric X-ray 
luminosity, =s 4 x 10 37 erg s -1 (M08). Note that most of the cascade 
power is emitted around ~10 GeV, i.e., within the band of Fermi 
GST. However, its sensitivit}0 allows to detect such a flare only for 
a > 30-hr duration. On the other hand, if there is a persistent emis- 
sion at a level a factor of 10 below that of the flare, it would be 
detected by Fermi GST within one year. 

In our model, we have assumed that the Compton losses of the 
cascade pairs dominate over the synchrotron losses. This requires 
that the magnetic field of the star is relatively weak, B < 5 G or 
so (see Fig.0 along the line of sight. The magnetic field of HDE 
226868 is unknown. In our model, it is required to be weaker than 
that of strongly-magn etized O stars, w here the surface dipole field 
can reach ~ 10 3 G (Wade et al. 2006). This is also pointed out by 
iBosch-Ramon et alj d2008h . We also assume that the pairs of the 
cascade quickly isotropise after their production, and produce 7- 
rays via Compton scattering locally, which will be the case even 
for a very weak field. 

An alternative explanation of the observed flare is that the 
primary TeV emission occurre d at a large distance from the 
black hole, e.g., along its jet (jPerucho & Bosch-Ramonl 1200 8t 
IBosch-Ramon & Khangulvanll2008| ), far enou gh to avoid the pair 
absor ption by the stellar photons (see also IBosch-Ramon et all 
2008). The observed TeV variability (A07) implies the emission 
region has a size < 10a. The distance fro m the black hole needs to 
be of the order of a few times a (fig. 2 in [Bednarek & Giovannellil 
|2007|) . Even if it is ~ la, it corresponds to ~ 2 x 10 6 /? g for M = 
1OM0. In the case of blazars, most of the y-ray emission is con- 
sider ed to come from d istances a few orders of magnitude smaller 
(e.g.. lSikoraet1d1ll994l) . On the other hand, this has been chal- 
lenged in some cases, especially for M87, where the TeV em ission 
has been argued to originate at > 120 pc dCheung et al.ll2007[) . For 
the black hole mass of 3 x 10 9 Mq, this corresponds to ~ 3 x 10 5 /? g , 
still substantially less than corresponding number for Cyg X-l. 

Note that some energy is transported to s uch distances in Cy g 
X-l in its jet, and it appears as radio emission (Stirling et al .1200 lh . 
Ho wever, the total radi o luminosity of Cyg X-l is only ~ 10 3 ' erg 
s 1 dFender et alj|2000l) . less than three orders of magnitude than 
the observed 0. 1-1 TeV power, and five orders of magnitude less 
then the intrinsic cascade power, see above. 



5 CONCLUSIONS 

We have studied opacity and acceleration models for the TeV flare 
observed from Cyg X-l by the MAGIC telescope. We have found 
that spatially extended pair cascades allow some TeV photons to 
escape pair absorption on the stellar photons. This can explain the 
observations provided the companion star in the system has a weak 
magnetic field. 



2 www-glast.slac.stanford.edu/software/IS/glast_lat_performance.htm 
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